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Abstract 

In this note, we relate the two well-known difficulties of Godunov schemes: the 
carbuncle phenomena in simulating high Mach number flow, and the inaccurate 
pressure profile in simulating low Mach number flow. We introduced two simple 
low- Mach- number modifications for the classical Roe flux to decrease the differ- 
ence between the acoustic and advection contributions of the numerical dissipation. 
While the first modification increases the local numerical dissipation, the second 
decreases it. The numerical tests on the double-Mach reflection problem show that 
both modifications eliminate the kinked Mach stem suffered by the original flux. 
These results suggest that, other than insufficient numerical dissipation near the 
shock front, the carbuncle phenomena is strongly relevant to the non-comparable 
acoustic and advection contributions of the numerical dissipation produced by Go- 
dunov schemes due to the low Mach number effect. 
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1 Introduction 



In the last decades, Godunov schemes are among the most successful meth- 
ods in simulating compressible flows involving shock waves and discontinuities 



181 ] . Godunov schemes utilize the solution of an Riemann solver at the face 



of computational cells as the numerical flux to introduce sufficient numeri- 
cal dissipation 5[. However, some low-dissipation Godunov schemes, such as 
those with Roe flux, may suffer from numerical instabilities at the strong shock 
front, known as the carbuncle phenomena, when simulating multi-dimensional 



high Mach number flows 15j. Although has not being fully understood, it is 
generally believed that this problem is caused by the insufficient numerical dis- 



sipation near the shock front 



15 
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8j. Therefore, up to now, almost 



all the cures proposed for Godunov schemes is trying to increase the numer- 
ical dissipation in the near shock-front region, or in the entire computational 



domain M, 
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|7j. 



In this short note, we propose to relate this problem to another well-known 
difficulty of Godunov schemes in the low Mach number limit 4j, y, |2| . The 
reason for this connection is that: for a shock wave in the solution of a multi- 
dimensional high Mach number flow, if it propagates in one direction of the 
Cartesian, or nearly Cartesian grid, the disturbances parallel to, especially 
behind, the front propagate in a low Mach number fashion. Actually, the 



alignment of shock front 



carbuncle phenomena 16, 



;o the grid is one of the typical scenarios of the 
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2 Roe flux and its low Mach number modifications 



For simplicity, we consider two-dimensional Euler compressible equation 

+ <9F(U) + 0G(U) = Q 

dt dx dy 



where U = (p, pu, pv, E) T , F(U) = [pu, pu 2 + p, puv, (E + p)u] T and G(U) = 
[pv, pvu, pv 2 + p, (E +p)v] T . This set of equations describes the conservation 
laws for mass density p, momentum density pv = (pu, pv) and total energy 
density E = pe + pv 2 /2, where e is the internal energy per unit mass. To 
close this set of equations, the ideal-gas equation of state p = (y — type with 
constant 7 is used. 

The classical Roe flux gives the following form of numerical flux 

Fj+i/2 = - (Fj+i + Fj) - -R^+i^lAj+i^lR^/^Uj+i - Uj), (2) 

where R and R _1 are the right and left eigenvector matrices of dF/dXJ and 
A the diagonal matrix formed with relevant eigenvalues: 

A lj2 = u, A 3>4 = u±c. (3) 

The first term on the right- hand- side of Eq. (j2J) is the central flux term, and the 
second term is the dissipative flux term. Since R and R _1 are only forward and 
backward coordinate transformations, the dissipative flux is merely dependent 
on A, and has two contributions: one is proportional to \u\, is called advection 
dissipation, the other proportional to \u ± c\, is called acoustic dissipation. 



The asymptotic analysis on Roe flux and general Godunov schemes 



,[3j show 



that, when M < 1, where M = \u\lc is the Mach number, or the acoustic 
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contribution of the numerical dissipation is much larger than that of the ad- 
vection, the dissipative flux term lead to pressure fluctuations of order O(M) 
at the cell face even if the initial data are well-prepared and contain only pres- 
sure disturbances of order 0(M 2 ). This phenomenon on inaccurate prediction 
of pressure profile is called low Mach number effect. One straightforward way 
to decrease the low Mach number effect is increasing the value of M in A 
to obtain comparable acoustic and advection contributions. This leads to two 
simple modifications of A to A' 

X' 12 — u, A3 4 = u ± min(0|w|, c), (4) 



and to A" 



A l,2 



u 



A 3,4 



u 



±c, 



u 



sgn(w) max 



(5) 



where is a positive number of order 0(1). The Roe flux with Eq. (T4|) is 



denoted as Roe-Mi and with Eq. (jSJ) as Roe-M2. Unlike Guillard and Viozat 



4] 



and Li et al. [12| . these modifications do not change the eigenvector matrix or 
the central flux term. Note that, though both modifications lead to comparable 
acoustic and advection contributions of numerical dissipation in the low Mach 
number region or direction of the flow, while the Roe-Mi flux decreases the 
local numerical dissipation, the Roe-M2 flux increases it. 



3 Simulation and discussion 

We test the calssical Roe, Roe-Mi and Roe-M2 fluxes, with 6 = 5, for a 

pi 

problem from Woodward and Colella [19( on the double Mach reflection of a 
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Mach 10 shock in air. The initial conditions are 



(p,u,v,p) 



(1.4,0,0,1) if y < 1.732 (x - 0.1667) 

(8,7.145,-4.125,116.8333) else 

and the final time is t = 0.2. The computational domain of this problem is 
[0,0] x [4,1]. Initially, the shock extends from the point x = 0.1667 at the 
bottom to the top of the computational domain. Along the bottom boundary, 
at y — 0, the region from x = to x = 0.1667 is always assigned post-shock 
conditions, whereas a reflecting wall condition is set from x = 0.1667 to x — 4. 
Inflow and outflow boundary conditions are applied at the left and right ends 
of the domain, respectively. The values at the top boundary are set to describe 
the exact motion of a Mach 10 shock. The calculations are carried out with 
the 5th-order WENO-Z [l] scheme, in which the Roe approximation is used 
for the characteristic decomposition at the cell faces and the 3rd-order TVD 
Runge-Kutta scheme is used for time integration 



11 



171 ] . Note that, in order to 



, |8j , no entropy fix is used in the 



exclude the influence of entropy condition 
computations. 



Figure [T] shows the density contours of the solution on a 480 x 120 grid at 
t = 0.2. Note that, a positivity-preserving scheme 6j is implemented to obtain 
numerically stable results for the classical Roe flux. It is observed that, contrast 
to that the classical Roe flux suffers from a kinked Mach stem, a typical 
configuration of the carbuncle phenomena, both Roe-Mi and Roe-M2 fluxes 
produce numerically stable and correct results. While Roe-M2 flux eliminates 
the kinked Mach stem by introducing extra advection dissipation parallel to 
the shock front, Roe-Mi flux achieves this by decreasing acoustic dissipation. 
Note that, since Roe-Mi decreases the overall local numerical dissipation near 
the shock front, it is unlikely that the kinked Mach stem is produced due to 
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insufficient numerical dissipation. Further numerical experiments show that, if 
the numerical dissipation is modified by only increasing the value of c in Eq. ([3]) 
i.e. only the acoustic dissipation increased, the simulation still suffers from the 
kinked Mach stem (not shown here). All these results clearly suggest that the 
carbuncle phenomena is strongly relevant to the non-comparable acoustic and 
advection contributions of the numerical dissipation produced by Godunov 
schemes due to the low Mach number effect. 
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Figure 1. Double-Mach reflection of a Mach 10 shock wave: 40 density contours from 
1.88783 to 20.9144 at t = 0.2, on a 480 x 120 grid with the (a) classical Roe, (b) 
Roe-Mi, and (c) Roe-M2 fluxes. 
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